Hypernetted-chain study of broken rotational symmetry states for the = 1/3 
fractional quantum Hall effect and other fractionally filled Landau levels 
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We investigate broken rotational symmetry (BRS) states for the fractional quantum Hall ef- 
fect (FQHE) at 1/3- filling of the valence Landau level (LL). Recent Monte Carlo calculations by 
Musaelian and Joynt [J. Phys.: Condens. Matter 8, L105 (1996)] suggest that Laughlin's state 
becomes unstable to a BRS state for some critical finite thickness value. We study in detail the 
properties of such state by performing a hypernetted-chain calculation that gives results in the ther- 
modynamic limit, complementing other methods which are limited to a finite number of particles. 
Our results indicate that while Laughlin's state is stable in the lowest LL, in higher LLs a BRS 
instability occurs, perhaps indicating the absence of FQHE at partial fillings of higher LLs. Possible 
connections to the newly discovered liquid crystalline phases in higher LLs are also discussed. 



PACS: 73.43.-f, 73.43.Lp, 73.43.Nq, 73.20.Mf, 64.70.Md. 



I. INTRODUCTION 

Recently, a plaethora of new phenomena has emerged 
in the transitional regions between different plateaus of 
the Hall conductanceoQ of Landau levels (LL) with in- 
dex L > 2. Near half- filling of the valence LL, ex- 
treme anisotropy has been measured in th& magneto- 
transport below temperatures ca. 100 mK,EJ"l3 accompa- 
nied by smooth non-linearities.El In addition, reentrant 
integer quantum Hall effect (RIQHE) regions with strik- 
ing breakdown features and ne,w phase transitions (pre- 
sumed to be quantum in|DriginEl) have been seen near 1/4 
filling of the valence LL.Q 

The anisotropic behavior has been attributed to the 
formation of a nematic phase of the two dimensional elec- 
tron system (2DES) which undergo-a. nematic to isotropic 
transition at higher temperatures-ElQ Similarly, the reen- 
trant regions are believed to be "bubble phases" similar 
to Wigner crystals, but with several electrons per bubble, 
or possibly new electronic hexatic states.lj 

The motivation of our work is to study these numer- 
ous liquid crystalline phases present in partially filled LLs 
by means of many-body trial wavefunctions with broken 
rotation symmetry (BRS), .similar to those proposed by 
Musaelian and Joynt (MJ)l3 in the context of the frac- 
tional quantum Hall effect (FQHE) (whichjare essentially 
generalizations of Laughlin's wavefunctiona) : 



*q(zi, ...,zn) 



(1) 



N 



Y[ (zi - Zj){zi - Zj ~ a){zi 



Zj +a) Z^i=i l^'l 



' i Hj is the complex 2D coordinate of j-th 



where Zj = 

electron and a is a complex number (we work in units of 
the magnetic length: Zq = h/eB = 1). This wavefunction 
represents a homogeneous liquid state with filling factor 
V — 1/3, lies entirely in the lowest LL (LLL), and for 
a ^ has nematic order (for a = we recover Laughlin's 
wavefunction, which is obviously isotropic). 



This wavefunction represents, therefore, a good start- 
ing point to consider nematic QH systems (even though 
the filling factor is "incorrect" for the newly discovered 
anisotropic states, see below for further analysis), by fa- 
cilitating the systematic study the energy dependence 
of BRS states for diverse physical parameters (LL in- 
dex, width of the 2DES, etc). In addition, the possi- 
bility of BRS states for v — 1/3 is intriguing by itselL 
since for a ^ Q highly damped low-energy modes exist ,E£I 
strongly modifying the dynamics and possibly suppress- 
ing the FQHE. 

In ReL || MJ investigated the possibility of BRS in a 
FQHE system using Eq. (|l|) and performing Monte Carlo 
(MC) simulations in a disc geometry. Their results sug- 
gest that the Laughlin fluid becomes unstable towards 
BRS states in quantum wells whose thickness exceeds a 
critical value depending on the electronic density. 

In this work we study the BRS state for 1 / 3-filling of 
the valence LL (i.e. ^ = M -I- 1/3 with M iirtegerljpy us- 
ing the hypernetted-chain (HNC) methodB&n This 
method allows us to compute physical quantities in the 
thermodynamic limit, without the limitations of using a 
finite number of particles that hinder other techniques, 
where the extrapolation to the thermodynamic limit is 
not totally unambiguous. We find that, contrary to MJ's 
results, the Laughlin state is stable in the LLL, whereas 
a BRS instability is possible in higher LLs (perhaps in- 
dicating, analogously to the arguments of Ref. |l^, why 
there is no ordinary FQHE at, e.g. v — 7/3). 

In Sec. H we present the basic theoretical calculations 
needed to determine the stability of an isotropic or BRS 
state. A detailed description of the HNC formalism in 
the context of the BRS wave function is given in Sec. 
HI. The results for the BRS state in the LLL and their 



extension in higher Landau levels are discussed in Sec. 
i^. Finally, in Sec. ^ we discuss our results, and analyze 
how they can be extended to more realistic filling factors. 
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II. BASIC THEORY 

In this work we propose to study the stability of differ- 
ent states by using trial wavefunctions like Eq. (|l|). We 
are interested, therefore, to calculate the energy in each 
of these states to find the optimum value for the sole free 
parameter a. 

We first consider the situation in the LLL, namely the 
state with i/ = 1/3. Since the BRS wave function is 
completely in the LLL the kinetic energy per particle is 
quenched at the lowest cyclotron energy 



1 (^ali^l^c.) _ 1. 



(2) 



where uJc = eB/m is the cyclotron frequency. The po- 
tential, or correlation, energy per electron is: 



En — 



N 



where V represents the electron-electron, electron- 
background, and background-background interaction; 
and g{r) is the (angle-dependent) pair distribution func- 
tion given by 

where r = r2 — ri. The following sum rule can be 
easilty proven p J d^r[g{r) — 1] = —1, and is a con- 
venient check for numerical procedures. For an ideal 
2D sample the interaction is a pure Coulomb poten- 
tial V{r) ~ e^/er , while in samples with finite thick- 
ness a reasonable choice is the Zhang Das Sarma (ZDS) 
potentialtS V{r) — e^/eV?"^ + A^, where A is of the order 
of the sample thickness. Alternatively, the correlation 
energy can be computed in reciprocal space: 
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(27r)2 



V{q) [5(q) - 1] , 



(5) 



where V{q) is the 2D Fourier transfor mQ (FT) of V{r) 
and S'(q) is the static structure factor: 



5(q)-l=pFT[.g(r)-l] 



(6) 



While both g(r) and S'(q) are angle-dependent (e.g. 
see Figs. and ||), because the interaction potential is 
centrally symmetric, the energy Ea depends only on the 
angle-averaged pair distribution function or static struc- 
ture factor defined as: 



de 

2^ 



5(r) , S{q) = 



27r 



5(q). (7) 



The determination of either the pair distribution func- 
tion or the structure factor is generally a complicated 
integral problem that needs to be solved for each LL. 
However, its is known that if transitions to other LLs 



are neglected (i.e. a single-LL approximation), g{r) and 
S{q) at higher LL are simply related to those at the LLL 
(L = 0) by means of a convolution or product respec- 
tively. We will apply this approximation (which, more- 
over, quenches the kinetic energy in higher LLs as well). 
It is then, sufficient to compute these distribution func- 
tions once in the LLL and then the correlation energy 
per electron is given by 



E^ = 



1 



(8) 



where Vcs{q) = V{q) [-^-l(9^/2)]^. Ll{z) are Laguerre 
polynomials, and S'(q) is calculated in the LLL [L — Q). 

In what follows we compute g{Y) and S'(q) in the LLL 
using the HNC method. 

III. THE HNC METHOD FOR THE BROKEN 
ROTATIONAL SYMMETRY STATE 

IntegmL equation techniques such as the HNC 
theoryOllj allow an accurate evaluation of the pair distri- 
bution function and related quantities associated with a 
Jastrow wave function. In particular they are extremely 
useful for calculations that are performed in the thermo- 
dynamic lipwt. They have been widely used in the study 
of classically and quantum fluids Jl3'll3 However the HNC 
method for the BRS wave function is a slightly different 
from that for the Laughlin wave function, since correla- 
tions and related quantities depend on both distance and 
relative angle between a pair of particles. 

The main quantity to be calculated in a HNC expan- 
sion is the pair distribution function g{v) [Eq. (^)], or 
equivalently the structure factor S'(q) [Eq. (0)]. These 
may then be used in conjunction with Eqs. (^, (|^), or 
(^) to determine the energy per electron for arbitrary 
values of the BRS parameter a or the 2D width A. 

Although the BRS wave function is a Fermi wave func- 
tion, its modulus square. 



|*a(2l, . . . ,2:Ar)| 
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(9) 



where u{z) = In |z|^-|-ln [z— a|^-|-ln [z-l-al^, can be viewed 
as a symmetric Jastrow wave function with pair correla- 
tions and single-particle terms. Theaefore it is possible 
to apply the Bose HNC formalismrJ In order to com- 
pute Eq. (^) one needs some small parameter in which 
to expand perturbatively (and re-sum a subset of dia- 
grams). For standard systems like Bose liquid "'He, the 
pair correlation is short-range and heals to 1 for large 
distances, therefore the function exp[u(ry)] — 1 provides 
a possible expansion parameter [note that in order to 
apply the Bose HNC expansion, the correlation (pseudo) 
potential has to satisfy the conditions: u{rij ^ 0) ^ —00 
and u{rij +00) — > 0]. In the case of the BRS wave 
function, the correlation (pseudo) potential is logarith- 
mically long-range, however it is possible to extend the 
method formally by splitting all quantities to compute 
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into a short- and long-range parts (see below). It can 
be shown that the pair distribution function can be ex- 
pressed as a series of cluster terms associated with linked 
diagrams and will be given from the following HNC equa- 
tions: 



X(ri2) = e''^'^i2)+N{r^2)+E{r^2) _ ^(^la) - 1 



(10) 



7V(ri2) = p I dVg X(ri3) • [X(r32) + A^(r32)] , (11) 



g(ri2) = 1 + X{ri2)+N{r 



12) 



(12) 



The quantities X{ri2) and N{ri2) represent the sum of 
the so-called composite and nodal diagrams respectively 
and E{ri2) is the sum of elementary diagrams. The 
generation of diagrams contributing to g(ri2) must go 
through a self-consistent procedure. As a first approx- 
imation (and a good one) we take the HNC/0 approxi- 
mation where the "0" denotes the neglect of elementary 
diagrams. The summation of the nodal diagrams iV(ri2) 
is easily performed in Fourier space. 

In order to handle the 2D logarithmic (pseudo) po- 
tential u(ri2), the standard procedure is to split it into 
short- and long-range parts: 



Zi(ri2) = Us(ri2) + Ui(i"i2) 



(13) 



with the nodal function iV(ri2) and the composite func- 
tion X{ri2) similarly split: 



N{ri2) = Nsiri2)-ui{ri2): 
X{ri2) = Xs{ri2) + u/(ri2) . 



(14) 
(15) 



This splitting is done subject to the following conditions: 

u(ri2) + iV(ri2) = u,(ri2) + Nsiru) , (16) 
7V(ri2) + X(ri2) = Nsiri2) + ^s(ri2) • (17) 

Given the particular form of (pseudo) potential for the 
BRS wave function, we choose to decompose u(ri2) into 
its short-range function (going to -oo for small distances 
and healing to for large distances) and its long-range 
counterpart in the following manner: 

U3(ri2) = -2A'o(Qri2) (18) 
-2Ko{Q\ri2~d\) - 2Xo(Q|ri2 + a\) , 

uj(ri2) = +2 [ln(ri2) + Ko{Qn2)] (19) 
+2[ln(|ri2 -a|) + Xo(Q|ri2 -a|)] 
+2 [ln(|ri2 + d\) + Ko{Q\ri2 + a\)] , 

where Ko{x) is the modified Bessel function, and Q isj-a, 
cut-off parameter of order 1. We recall that the 2D FTxB 
of u/(ri2) is: 



w/(q) 



47rQ2 



q^iq^ + Q2) 



1 



(20) 



The final set of equations is solved by initially set- 
ting A^s(ri2) = in Eq. (p^), then obtaining Xs{q) = 



FT[Xs(ri2)] which can be used to compute X(q) = 
^s(q) + ""((l)- Using the convolution theorem we find 
iV(q) = pX(q)2/[l - pX{q)] and easily obtain Ns{q) = 
N{q) + ui{ci). The last step is to perform an inverse 2D 
FT on Ns{l) to obtain the new A^s(ri2). This procedure 
is repeated until a desired accuracy is reached. After 
convergence the pair distribution function is given by 

g(ri2) = l + X,(ri2)+iV,(ri2). (21) 

Simultaneously, the static structure factor is given by 



S{q) 



1 



l,(q)+iV,(q) 



(22) 



The computation of such functions allows us to find the 
interaction energy and other related quantities. 

IV. RESULTS AND DISCUSSIONS 

In the present work we applied the HNC theory to 
study the BRS state at filling 1/3 of an arbitrary LL (in 
the single-LL approximation). For the sake of simplicity 
we neglected the elementary diagrams (i.e. the so-called 
HNC/0 approximation). This allows us to determine to 
a reasonable accuracy the pair distribution function and 
the static structure factor. In order to compare the a = 
(Laughlin) state with the a 7^ (BRS) state we studied 
the properties of the BRS wave function for several as 
with magnitudes between to 3 (in general a — \a\ e^^" , 
without loosing generality we considered only 9a =0). 

A. Pair distribution function and structure factor 

In Fig. |l| we plot the pair distribution function f;(r) 
for a = 2 (top, center panels), and the angle-averaged 
pair distribution function g{r) corresponding to a = 0, 
1, 2 and 3 (bottom panel). It is interesting to note the 
splitting of the triple node at the origin, the noticeable 
angle-dependence of ^(r), and the change in the small-r 
behavior of g{r) which switches from oc (for a = 0) 
to oc as a is increased. In addition, we mention the 
following generic properties more or less valid for any of 
the a-s we have considered: 

1. For a 7^ there is at least one pair distribution 
function that has an additional zero (besides the 
zero at r = 0) at inter-particle distance r — a and 
angle 6 = 6a,0a + tt {9a = in this case); 

2. For a 7^ there are special inter-particle distances, 
r (besides the zero at the origin) where all pair dis- 
tribution functions cross, irrespective of their angle 
9 dependence. 

3. Extremely interesting is the behavior of the angle- 
averaged pair distribution function 'g{r) as a func- 
tion of a. One notes that the major peak of ^(r) 
simply shifts to larger distances (without any size- 
able change in its height) as a is increased. For 
smaller distances g{r) starts to develop a shoulder 



3 



that is quite visible for a = 3 contrary to what 
seen in Ref. || for a slightly larger a = 3.2, where 
the shoulder should had been even larger. 

4. For small-r, ^(r) has almost no angular depen- 
dence, and for a ^ 0, g{r « 0,6*) ~ Ca for 
< r < 0.5, where Ca ^ 0.026 a^ ^ when a = 0, 
g{r « 0) oc as expected. These results derive 
immediately from the 1-fold vanishing of the BRS 
wave function when two electrons come close, as 
opposed to the 3-fold vanishing of Laughlin's wave 
function, since for small distances only two-body 
correlations are important. The absence of angu- 
lar dependence on g(r, 9) for small r can be eas- 
ily understood by noting that in the small-r limit: 
gir w 0,6') cx exp[ws(r « 0,6^)], where Us{r,6) is 
given from Eq. and does not have any angular 
dependence. By recalling that limr-,o Ko{Qr) — 
-ln(Qr/2) - 7, where 7 = 0.5772... is the Euler's 
constant one can easily understand why g{r,9) has 
a quadratic dependence on r and not any angu- 
lar dependence for small r and values of a ^ 0. 
Although such quadratic dependence at short r is 
also characteristic for a Wigner crystal state at such 
filling factoiO we note that the BRS state does not 
represent a crystalline state and the overall pair 
distribution function of the BRS state is strikingly 
different from the pair distribution function of the 
Wigner crystal state. 

In Fig. H we plot the static structure factor 5(q) for 
a — 2 (top, center panels), and the angle-averaged static 
structure factor S{q) corresponding to a = 0, 1,2 and 3 
(bottom panel) . The most important feature is the emer- 
gence of peaks in S{q) characteristic of a nematic struc- 
ture. Broadly speaking, the major peak of the S{q) shifts 
to smaller q and its height raises when a is increased, with 
no significant change in the small-g behavior. 

B. Energy of BRS states 

One can compute the correlation energy per particle 
either directly from Eqs. (||), (||), or (||) to determine 
the energy per electron for arbitrary values of the BRS 
parameter a, the 2D system width A, and Landau level 
index L. The following simplified formula can be used in 
view of Eq. (|^): 

^ r°° ~ — 
EtW = ^l dqqV{q,X)[LUY)]'iS{q)-l], (23) 

where V{q,X) = {2TTe'^/eq) e^{-Xq) is the 2D FT of 
the ZDS interaction potential£3 In addition to allowing 
strai ght forward calculations to be extended to any LL, 
Eq. ( p3| ) permits higher numerical accuracy on the cal- 
culation of Ea since S{q) saturates exponentially to 1 for 
relatively small values of q as compared to g{r). 

Figure || shows the energy difference between BRS 
states with a = 1, 2 and 3, and the isotropic state with 



a = 0. Our findings indicate that in the LLL {L = 0) 
the Laughlin state is stable for any A, since all a ^ 
states have higher energies (top panel), contrary to prior 
resultaa that the BRS state for a = 1 has a lower en- 
ergy than the Laughlin state if one considers the ZDS 
potential, with A > Ac = 4.1 ± 1.5. 

The situation changes drammatically in higher LLs 
{L > 1). For small A BRS states have lower energies 
and the incompressible Laughlin-like state is unstable 
(see lower panels of Fig. ||) . The presence of this insta- 
bility towards a BRS state may be related to the absence 
of FQHE states in higher LLs, since for a ^ highly 
damped_Low-energy modes exist in the resulting nematic 
system.EJ It is worth noting that for A ^ 1 the high- 
est investigated a has the lowest energy. In this regime, 
we are therefore unable to determine the optimal state 
(even within this familty of trial wavefunctions). This 
BRS instability may be indicative of a transition towards 
a completely different state (e.g. as in Ref. |l^). 

At this point it is important to comment on how pre- 
cise our determination of these energy differences is. The 
reader should note that the HNC /G approximation is es- 
sentially a variational method that always gives an en- 
ergy that constitutes an upper bound to the exact ground 
state energy.Eil For example: for the Laughlin state with 

= 1/3 and A = 0, HNC/O gives an adimensional cor- 
relation energy of -0.40&- while the "exact" value (deter- 
mined by Monte Carloc2l) is -0.410. Similar errors (ca. 
1%) will be present for a ^ as well. While an error of 
this magnitude seems to be of the same order as, or big- 
ger than, the energy differences we are interested in, we 
remark that these are not uncorrelated errors but system- 
atic deviations due to the nature of the approximations 
used, and energy differences are considerably more pr|e=| 
cise. Preliminary results using Monte Carlo simulationsc3 
for a handful of cases indicate that energy differences are, 
indeed, significant. 

V. CONCLUSIONS, EXTENSIONS AND 
FURTHER DEVELOPMENTS 

In conclusion, we applied the HNC theory to study 
possible BRS states in a 1/3-filled LL. We find that the 
isotropic Laughlin state is stable in the LLL for real- 
istic interaction potentials. In higher LLs, instabilities 
towardg-a BRS state are possible. Since BRS states are 
gaplessE3 this may be a simple explanation why no FQHE 
was observed for 1 /3-filled higher LLs. One caveat is that 
the magnitude of the energy differences obtained is com- 
parable to the absolute accuracy in the determination of 
individual energies by the HNC method. Although we 
believe that energy differences are yielded more precisely 
than the energies themselves, these results need confirma- 
tion by alternative (albeit more time-consuming) meth- 
ods. Monte Carlo simulations with la*ge number of elec- 
trons are currently being performed.E3 

While these results are by themselves compelling, 
the connection to recent obsevations of liquid cristalline 
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phases in half- and quarter-filled LLs requires more so- 
phisticated methods. One possibility is to generalize 
MJ's approach to composite fermion (CF) states,E3 e.g. 

^]I\z^,...,zn)= (24) 

Pl Wizj - Zk + a){zj -Zk-a) Det [0k(rO]fc<fc^ , 

]<k 

where 0k (fi) are plane waves of CFs, filled for k<kF — 
(47rp)^/^, and projects onto the L"^ LL. This wave- 
function is an obvious starting point to study the ne- 
matic quantum Hall liquid crystals at half filling. For 
the RIQHE observed near 1/4 filling, similar generaliza- 
tions are possible. The presence of the Slater determinant 
in the trial (in addition to the Jastrow factors) implies 
the ipfrfn to use the considerably more complex Fermi 
HNQlilo. Calculations are under way for trial states of 
these form.E^ 
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We use the standard convention for the 2D FT: /(q) = 
Jd\ expHq-r] /(r), /(r) = / d^q/i27vf exp[jq.r] /(q). 
For half-filled LLs, more appropriate trial wavefua|Ctions 
may start with a composite fermion type statecJ with 
anisotropic Jastrow factors [similar to Eq. (jl|)] and an 
isotropic Slater determinant. These type of wavefunctions 
requirie-ijthe use of the considerably more complex Fermi 
HNC,yy see comments in Sec. 0. 
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FIGURES 
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FIG. 1. Pair distribution function for the BRS state at 
v = 1/3. (a) a = 2, surface plot of g{r,9) (the surface for 
y < was removed for clarity); (b) a = 2, dotted lines: 
g{r,6) for various G [0, 27r], full line: angle averaged g{r); 
(c) Angle averaged g{r) for a = 0, 1, 2 and 3 (0 and 1 are vir- 
tually identical). Note the discrete nodes of g{r, 6) at r = a, 
6 = 6a,0a + IT {6a = in this case). Calculations were per- 
formed in the HNC/0 approximation. 
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FIG. 2. Static structure factor for the BRS state at 
v = 1/3. (a) a = 2, surface plot of S{r,6) (the surface for 
% < was removed for clarity); (b) a = 2, dotted lines: 
S{q, 9q) for various 9q £ [0, 27r], full line: angle averaged S{q); 
(c) Angle averaged S{q) for a = 0, 1, 2 and 3 (0 and 1 are 
virtually identical). Note the presence of peaks in S'(q) con- 
sistent with a nematic structure. Calculations were performed 
in the HNC/0 approximation. 
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FIG. 3. Energy per particle in BRS states with a — 1,2 and 3 relative to the isotropic (a = 0) state: AEa{X) — Ea{\) — Eo{X) 
for various Landau levels L as functions of the short distance cut-off A [Eq. (|2^)]. Energies are in units of e^/(eZo). Note that in 
the LLL (L = 0) BRS states are always higher in energy, whereas in higher LLs {L > 1) there are ranges of A for which BRS 
states are favorable. 
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